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Abstract 



We extend the McLerran-Venugopalan model for the gluon distribution func- 
tions of very large nuclei to larger values of the longitudinal momentum frac- 
tion Xp. Because gluons with larger values of Xp begin to resolve the lon- 
gitudinal structure of the nucleus, we find that it is necessary to set up a 
fully three-dimensional formalism for performing the calculation. We obtain 
a relatively compact expression for the gluon number density provided that 
the nucleus is sufficiently large and consists of color-neutral nucleons. Our ex- 
pressions for the gluon number density saturate at small transverse momenta. 
The nuclear dependence we obtain is such that the number of gluons increases 
more slowly than the number of nucleons is increased. 

24.85. +p, 12.38.Bx 
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I. INTRODUCTION 



The recent construction and commissioning of the Brookhaven Relativistic Heavy Ion 
CoUider has led to a renewed interest in the properties of heavy nuclei in recent years. A 
considerable amount of fruitful work has been done on classical and semiclassical descrip- 
tions of the physics involved |l|-|T9|. In particular, the McLerran-Venugopalan (MV) model 
provides a framework for calculating the gluon distribution functions for very large 
nuclei at very small values of the longitudinal momentum fraction x^. What McLerran and 
Venugopalan realized is that at sufficiently small Xp, the gluons are unable to resolve the 
longitudinal structure of the nucleus, meaning that many quarks contribute to the color field 
at each value of the (transverse) position x. This large charge per unit area provides 
the scale at which the strong coupling is evaluated [Q]. Thus, if ^ ^qcd; ^ classical 
treatment ought to provide a reasonable description. 

Recently, we pointed out that the infrared divergences which appear in the MV model 
may be cured by incorporating the effects of confinement 0. That is, we observe that 
nucleons display no net color charge: individual quarks are confined inside the nucleons, 
whose radius is a ~ Aqqj-,. As a consequence, we expect that there should not be long range 
(S> a) correlations between quarks. Strong correlations between quarks occur only when 
we probe at short distance scales. These considerations may be phrased as a mathematical 
constraint on the form of the two-point charge density correlation function. 

The calculations presented in Refs. ||T]-^ all assume Xp to be small enough so that the 
gluons do not probe the longitudinal structure of the Lorentz-contracted nucleus which they 
see. Effectively, then, the relevant geometry is two-dimensional, with the source exactly on 
the light cone. In this work we extend the MV model to larger values of the longitudinal 
momentum fraction Xp. In this regime, the gluons begin to resolve the longitudinal structure 
of the nucleus: therefore, we develop a fully three-dimensional framework using a source that 
is not quite aligned with the light cone. In order to deal with the complications which arise 
as a result, we must rely heavily on the fact that we consider a large nucleus of radius a 
which consists of color-neutral nucleons. 

The remainder of this paper is organized as follows. In Sec. y we present our conventions 
for writing down the classical Yang- Mills equations and describe the (slightly) off-light-cone 
source which will be the foundation of our calculation. We show that the natural (order 
unity) variables to describe the nucleus are essentially those in the nuclear rest frame, even 
in the limit (3 —>■ 1. In Sec. Ill we set up the framework for determining the gluon number 
density in the Weizsacker- Williams approximation. In this section we introduce the two- 
point charge density correlation function, and review the color-neutrality condition which 
it must satisfy. Sec. |rV| contains a discussion of the requirements which must be satisfied in 
order for our approximations to be valid. The meat of our calculation is contained in Sec. 0, 
where we begin with the solution for the vector potential in the covariant gauge, perform the 
transformation to light-cone gauge, and determine the gluon number density. Additional 
details of this calculation are found in Appendix ^. We illustrate our results for the gluon 
number density with the help of a power-law model for the correlation function in Sec. [VI| . 
The integrals which arise in connection with this model are presented in Appendix 0. Finally, 
Sec. IVIII contains our conclusions. 
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II. THE CLASSICAL YANG-MILLS EQUATIONS AND SOURCE 



In this section we present the conventions which we use in writing down the Yang-Mills 
equations and source for a nucleus which moves down the z axis with a speed /3 < 1. The 
motivation behind the choices we have made is to ensure that all "unknown" quantities are 
of order unity, with all powers of the small and large parameters explicitly written out. This 
will make the approximations which we will have to make later on more transparent. It will 
also make the /? — > 1 limit which we take at the end obvious. 

We begin with the classical Yang-Mills equations, which we write as 

D^F^^ = gr, (2.1) 

where we have employed matrix form, i.e. J'^ = T'^J"'", etc. The T*^ are the normalized 
Hermitian generators of SU(A'"c) in the fundamental representation, satisfying 2Tr(T"T^) = 
^ab^ The covariant derivative is 

D^F^" = d^F^" - ig [A^, F^''] (2.2) 

and the field strength reads 

F^"" = d^A" - a^A^ - ig [a^", A"] . (2.3) 

The conventions contained in Eqs. (|2.1| )- (|2.3| ) ensure that all powers of the strong coupling 
constant g are explicit, with no hidden (yf-dependence. 

We now turn to the source appearing in ( p.l|) , starting in the rest frame of the nucleus. 
Using the subscript "r" to denote rest frame quantities, the current takes on the simple form 

jO = p{-z,, Xr)] Jl = = = 0. (2.4) 

The color charge density p = T"'p"' is a spherically symmetric function which is non-zero 
over a region of size R, the radius of the nucleus. Since in the lab frame we want the nucleus 
to be moving along the +z axis, it is convenient to use —Zr for the longitudinal coordinate 
in Eq. ( p.4| ). The transverse coordinates Xr and yr form a two- vector which we write in 
bold-face: Xj.. In terms of the light-cone coordinates^ = —X:^ = (a;° ± x^)/v^, Eq. ( p.4|) 
may be written as 

-^r = J; = 7lP(7l(^;-^.^).^0' = 0- (2-5) 

The net color charge of the nucleus is zero: 

J dZr<fXr p{ — Zr,Xr) = 0. (2.6) 



^ Our metric has the signature (—,+,+,+). Thus, the scalar product in Ught-cone coordinates 
reads g^x^ = —q~^x~ — q~x~^ + q ■ x. We will think of as the time, and x~ as the longitudinal 
distance. 
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The nucleus is not a homogeneous sphere of color charge: it has substructure. Because 
of confinement, there are smaller regions of size a ~ Aq^j-, within the volume occupied by 
the nucleus for which the total color charge also vanishes. These regions correspond to the 
nucleons. For a large nucleus, a/i? ~ A~^^^ -C 1. 

Boosting to the lab frame, where the nucleus moves along the +z axis with a speed (3, 
Eq. (|2.5|) becomes 



e 



J^ = -p[lx--lx+,x); J' = 2'^^-^ J = 0^ (2-7) 



where we have defined 



1 + P 



(2.8) 



Viewed in the lab frame, the nucleus is Lorentz-contracted to a thickness of order -R/7, 
where, as usual, 7 = (1—/?^)"^/^. For large boosts the charge density function is non-zero 
only when x~ ~ eR. Hence, the longitudinal argument in Eq. ( |2.7D is really of order R, 
leading us to define the new longitudinal variable 

1 _ e , 

X|| = - X X . (2.9) 

The current appearing in Eq. ( p.7|) is a function of ( The use of the notation "x" 

is suggestive of the fact that, in terms of the natural (order unity) variables, the functions 
describing the nucleus are still spherical! In fact, Xy is just the (unboosted) longitudinal 
coordinate from the rest frame. The advantage of using Xy instead of x~ should be obvious: 
to take the (3 1 {e ^ 0) limit for quantities written in terms of Xy is trivial, whereas if 
the same quantities were written in terms of x~ instead, we would have to be careful to hold 
x~ /e fixed. 

Because the choice made in (E.4|), we define 



g ■ X = -g|[X|[ + qr • a;. (2.10) 

We will also use the notation d^x = dxy (Px. All of the functions we will be dealing with 
will depend only on the particular combination of x"*" and x~ appearing in Eq. (T^). Thus, 
we have the replacements 

The divergence of the current in Eq. ( |2.7| ) vanishes: d^J^ = 0. In QCD, however, we 
require that the current be covariantly conserved: D^J^ = 0. Since we are at weak coupling, 
we may work iteratively. That is, we first solve the Yang-Mills equations using ( p^.7| ) for 
the source. The resulting solution will violate the covariant conservation condition by an 
amount of order g"^. The value of Z^^J'^ which is obtained could, in principle, be used 
to correct the current to this order given the non-Abelian equivalent of the Lorentz force 
equation to provide information on how the color of the quarks making up the source changes 
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upon emission of a gluon. Since we are assuming that g"^ <^ 1, we will simply drop these 
contributions. 

As suggested by Eq. ( p.lQ ), we will use the notation for the component of momentum 
conjugate to To make the connection with the longitudinal momentum fraction x^, let 
the gluon carry a momentum and the nucleon a momentum . Then, 

g+ £g+ 

where m is the nucleon mass. However, 

q+ ^id_^-d«^^. (2.13) 

e e 

Thus, we conclude that 

= ^. (2.14) 

m 

Once more note the advantage of the rest-frame variables over the light-cone variables: in 
terms of q~^, we would have to take the e — limit with the caveat that the combination 
eq~^ is held fixed. No such complication arises when we use gy instead. 



III. COUNTING GLUONS 



expression reads |]20 
dN g^ 



Next, we turn to the formula for the gluon number density. Recall that the standard 
3 in 

r dx- r dx'- I d^x f d^x' e-'''^^''~-'''~^e"'<'^-'''\A'}{x-,x)A1{x'-,x')). 
dq d Q 47r J—oo J—oc J J 

(3.1) 

Eq. ( p.l[ ) is written in terms of the light-cone variables x~ and g"*", and the light-cone gauge 
vector potential Ai. The light-cone gauge has favored status with respect to the intuitive 



picture of the parton model p0|-p3|: thus we continue to use the light-cone gauge even for 
a source which moves at less than the speed of light. Based on the discussion of Sec. ||, 
however, we wish to employ the "new" longitudinal variables and gy; therefore, we write 



dN 



d^x I e'^''-(^-^')(v4^(x)4"(f' )) 



= ^ jd'xj d'x'e'''<^-^'\Atiq,;x)A^i-q,;x')), 

where 

/oo 
dx^^ e"*^"""'" A{x). 
-oo 

In the limit e —>■ 0, Eq. (p. 21) reduces to the previous result, Eq. (|3.1|). 



(3.2) 



(3.3) 
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Eq. ( |3.2| ) produces a gluon number density which is differential not only in gy, but in 
the transverse momentum as well. To obtain the usual gluon structure function resolved at 
the scale , we simply supply the trivial factor of m required to convert g,, into Xp and 
integrate ( p.2|) over all transverse momenta less than or equal to Q: 

gA{Xp,Q^)=f d^q-^^. (3.4) 
J\q\<Q dXpd^q 

Our classical approximation to the quantum average represented by the angled brackets 
appearing on the right-hand-side of Eq. ( |3.2| ) consists of performing an ensemble average 
with a Gaussian weight. We parameterize the two-point charge density correlation function 
by 

{p%x)p''{x')) = 6"' K^s(^^^^^V{x - x'). (3.5) 

The functions S and V appearing in this definition encode two different aspects of the physics 
of the nucleus. The color-neutrality condition developed in Ref. ^ imposes the following 
constraint on V: 



d'^xV{x) = 0. (3.6) 

In terms of the Fourier-transformed function V, this constraint reads 

:D(0,0) = 0. (3.7) 

When the color-neutrality condition is satisfied, the function V contains an intrinsic scale, 
reflecting the minimum size of the region for which Eq. ( ^.6] ) is approximately true. Because 
of confinement, we expect this scale to be roughly the nucleon radius a ~ Aqqj-,, reflecting the 
fact that points within different nucleons ought to be (largely) uncorrelated. On the other 
hand, the function S(^^{x+x')j , which depends on the center-of-mass coordinate, should be 

non- negligible over a region of size R ~ A^^^a, the radius of the entire nucleus. We choose 
to normalize S so that its total integral simply gives the volume of the nucleus: 

d^^ 5(S ) = V. (3.8) 

We have written V in Eq. ( p. 81) rather than the |vri?^ pertaining to a spherical nucleus to 
maintain generality and to aid in making connection to Refs. []5|,^. The detailed forms of 
the functions S and V depend upon aspects of non-perturbative QCD which are poorly 
understood. However, the requirements specified above account for the relevant physics: 
we have a correlation function which takes on non-trivial values only for points which are 
close enough together to lie inside a single nucleon {\x — x'\ ^ a) and which are centered 
anywhere inside the nucleus (||(x + x')| ^ R). The interplay between S and V will be 
crucial in helping us to organize our calculation in powers of a/R. 

In addition to satisfying the color neutrality condition, T>{x — x') should contain a term 
which goes like S^{x — x') [p^-^pAj. The presence of such a contribution leads to a point-like 
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behaviour in the gluon number density at large 0, consistent with the physics 
of asymptotic freedom. This term arises from the self-correlation of the quarks. Quite 
generally, then, we expect the form of T>{x — x') to be 



V{x-x')=6'{x-x')-C{x-x'), (3.9) 

where C{x — x') is a reasonably smooth function parameterizing the mutual correlations 
between pairs of quarks. Because C{x—x') describes the structure of a color-neutral nucleon, 
it must have unit integral and possess non-trivial values only when |x — af' | ^ a ~ Aqqj-,. In 
momentum space, Eq. ( |3.9D reads 

V{q) = l-Ciq). (3.10) 



The color neutrality condition ( |3.7| ) implies that C(0 ) = 1. Furthermore, for asymptotically 
large q, C{q) —>■ 0, since in position space C is reasonably smooth and has a finite region of 
support. 

Finally, the only quantity appearing in Eq. (p.5|) yet to be specified is K^. We determine 
by integrating the trace of ( p.5| ) and replacing V^A) by 6^{A): the result should be 
3ACf for a nucleus containing A nucleons. Thus 

,_3AC^ 1 _3A 1 



IV. REGION OF VALIDITY 

All of the machinery assembled in the previous two sections has been geared towards 
performing a classical computation of the vector potential associated with a color charge 
moving down the z axis with a speed P near, but not equal to, the speed of light. The vector 
potential is then translated into a gluon number density in the spirit of the Weizsacker- 
Williams approximation. In this subsection we will consider the conditions which must be 
satisfied in order for this treatment to be valid. 

Firstly, we need the coupling as = g'^/4:7i to be weak. When <^ 1, we have the pos- 
sibility that the quantum corrections will be small, making the classical result a reasonable 
approximation to the full result. Several years ago, McLerran and Venugopalan |^ observed 
that for a very large nucleus or at very small values of x^ the density of quarks and gluons 
per unit area per unit rapidity is large. When this density is much larger than Aqqj^, we 
expect as to be weak 

The large density of color charge facilitates the classical treatment in a second fashion: 
when there a large number of charges contributing to the charge density at each point, the 
total will (typically) be in a large representation of the gauge group. Thus, we may treat 
the source classically. In addition, the large number of quarks justifies the use of a Gaussian 
weight for the ensemble average via the central limit theorem. 

Under what conditions do we see a high color charge density? And when do a large 
number of charges contribute? The answers to these questions depend upon the scale at 
which the gluon distribution is being probed (see Fig. |I|). In the original MV treatment 
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|T|-P, it was assumed that was "small." In this case, "small" means that the longitudinal 
scale resolved by the gluons is larger than the Lorentz-contracted thickness of the nucleus 
which they see: 

1 

— >R, or x^< . (4.1) 

gii ma 

All of the quarks at a given transverse position x contribute to generating the gluon field 
measured at the values of x^ indicated in Eq. ( [4.1|) . This leads to a color charge per unit 
area of 

Asking that this density be ^ Aq^^j^ ~ so that the coupling be weak leads to the 
condition 

> JL_. (4.3) 

Refs. ||^-|^ restrict the allowable to the region 

Aqcd ^q'^ SA^/^CM^cD- (4.4) 



The lower limit in ( [4.4|) comes from the requirement that the gluons probe distances small 
compared to the nucleon radius. On the other hand, if is too large, an insufficient number 
of quarks will contribute, no matter how big the color charge density is. Since the amount of 
transverse area probed by a gluon with transverse momentum q^ is about vr/q^, we conclude 
that for the charge per unit area given in ( ^.2| ), the gluon sees an amount of charge equal 
to SA^/^Ci^AQQ^/qr^. Asking that this be much greater than one leads to the upper limit 
in ( [4.41 ). The requirements of Eqs. ( [4.1| ) and ( [4.4| ) restrict the range of validity of the original 
MV model to the region labelled "A" on Fig. ||. 

Because the strong coupling is evaluated at rather than q^, it ought to be possible 
to relax the lower limit in ( [4.4| ), provided that the theory is infrared finite. All that is 
required is a framework which captures the key consequence of confinement, namely the 
fact that when viewed on large (^ a) distance scales, the nucleons are color neutral. This 
observation leads to the color neutrality condition ( p.6|) to be imposed on the two-point 
charge density correlation function 0. Not only does the color neutrality condition make 
the theory infrared finite, but it also limits the amount of color charge being probed as 
g ^ 0: beyond about q^ ~ 1/a^ the net charge drops as complete color-neutral nucleons 
are probed. The net effect is to reduce /t^ by the factor (ag)^ from the value given in ( [4.2|) , 
leading to the less-stringent lower limit 



q'^-^^lcT.- (4.5) 



The region labelled "B" in Fig. |l| represents the additional range of validity obtained in 
infrared finite theories by replacing the lower limit of (|4.4|) with 
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Now we turn to the main goal of this paper, the relaxation of the condition (|4.1| ) on 
Xp. At larger values of x^, the gluons are able to probe shorter longitudinal distance scales: 
they no longer see the entire thickness of the Lorentz-contracted nucleus. The fraction of 
this thickness which they do see is roughly l/q\\R ~ l/mXpA^^^a. Hence, the value of k"^ 
obtained in Eq. ( [4 .21 ) is reduced by this factor, and the coupling is weak only for 

X, < (4.6) 
vrma 

Taking ma ~ 5, Eq. (|4.6| ) implies an upper limit of Xp ^ 0.25. At larger values of x^, too 
little of the nucleus is seen by the gluons in order for the charge density to be large and the 
coupling weak, independent of how large we imagine the nucleus to be. 

Eq. (O) is not the final word, however. As noted above, when becomes small, the 



effective charge density is reduced since we begin to see color neutral nucleons. Thus, we 
should further reduce k"^ by the factor (aq)'^ in this region, leading to the lower limit 



q^ > nma 



A^CD, (4.7) 



valid whenever x^ is "large" {i.e. when x^ is larger than the value given in Eq. ([4.1|) ). 
Likewise, when q^ is made too large, not enough charge is probed. The upper limit is 
reduced from the value given in Eq. ( [4.4|) by a factor of l/q\\R, to 

x,g^<^A|cD- (4.8) 
ma 

Taken together, the constraints ( [i.6[ ), ( [4.7D , and ( [4. 81 ) allow us to extend the computation 
of the gluon number density into the region labelled "C" on Fig. |l[ 

Finally, we note that the eikonal approximation which we are using also tells us that q^ 
and Xp cannot get too large: that is, we are ignoring nuclear recoil effects. 



V. GLUON NUMBER DENSITY 

We now turn to the computation of the gluon number density within the 3-dimensional 
framework described in Sees. || and |T|. In order to obtain our result, we will have to rely on 
both color-neutrality and the large nucleus approximation extensively. Our final expression 
reduces to the MV result of Refs. in the limit x^, — > 0, but only if the nucleus is assumed 
to have cylindrical geometry. 

Our calculation has three stages: first, we obtain the solution for the vector potential 
in the covariant gauge. Next, we transform that solution to the light-cone gauge. Finally, 
we use the light-cone gauge solution along with the correlation function ( p.5|) to obtain the 
gluon number density from Eq. (|3.2| ). 



A. Covariant Gauge Vector Potential 



The most efficient route in performing our calculation begins by imagining the situation 
in the rest frame of the nucleus, where we consider a static distribution of color charge, 
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Eq. (|2.4|) . In this frame we have the "obvious" time-independent Coulomb solution for the 
vector potential. Since only ^ 0, we have OqA^ = d ■ A = 0, that is, the Coulomb solution 
is the same as the covariant gauge solution. When we boost to the lab frame then, it is 
natural to begin with the covariant gauge solution. 

The Yang-Mills equations in the covariant gauge readFl 



2d+d^)A'' = gr + 2ig A^ , d^A'' 



- W 



A 



A", A" 



(5.1) 



In order to deal with these equations, we must assume not only that the source has the 
form indicated in Eq. ( ^.7| ), but also that we are in the weak-coupling regime, -C 1. In 
particular, we assume that the commutator terms appearing in (|5.1|) are negligible, leaving 
the simpler equations 

{V' + d^)A''{x)=gr{x). (5.2) 

The operator appearing in Eq. ( ^.2|) is simply the Laplacian in 3-dimensions. These equations 
are solved in the usual manner by introducing the Greens function G{x) which satisfies the 
equation 

iV^ + dl)Gix)=5%x). 



Passing to momentum space, we find that 

Giq) = 



-1 



The Fourier transform used to obtain ( |5.4| ) is easily inverted, producing 



(5.3) 



(5.4) 



(5.5) 



XT, 



Actually, because of the unequal treatment of the longitudinal and transverse variables when 
we transform to the light-cone gauge, the following mixed representation. 



Gix 



(5.6) 



47r2 2q 

which is obtained by inverting only the longitudinal part of the transform, will prove to be 
especially useful. Note that in Eq. ( p.6| ), as elsewhere in this paper, q means \q\. 

Independent of how we choose to write down the Greens function, the solution to 
with the source ( p.7| ) reads 

A^{x ) = - g J d^x' G{x — x')p{x') 



A {x\ 
A^{x\ 



e'A-^ix 



i-2 
2 

0. 



(5.7) 



^From this point forth, we will use a tilde to distinguish the vector potential in the covariant 
gauge from the vector potential in the light-cone gauge. 
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B. Transformation to the Light Cone Gauge 



At this stage, we are ready to perform the transformation to the hght-cone gauge. In a 
non-Abehan theory, we may parameterize the gauge transformation as 



(5.8) 



Since the potentials we are deahng with are functions of x only, we expect that U will 
also depend on x only. Thus, the requirement that the new gauge be the light-cone gauge 
becomes 



a||U(x) = ige'[J{x)A+{x), 



(5.9) 



where we have replaced d- by in accordance with (|2.11|) . The solution to Eq. ( |5.9| ) is the 
path-ordered exponential 



U(x ) = Pexp 



(5.10) 



In Eq. (|5.10| ) we have introduced the shorthand notation 

J —OQ J —OQ J —OO J —OQ 



dy\\ 



(5.11) 



to indicate the ordered integration region Xy >y\\i > y\\2 > ■ ■ ■ > > — oo. Eq. ( |5.7| ) tells 
us that A~^ is naturally of order 1/e: therefore, all of the terms in the sum on the right hand 
side of (|5.1CI|) are of order unity. Introducing the expression for A'^ into the expression for 
U(x ) produces 



U(x) = "Pexp 



(5.12) 



Inserting this into Eq. (|5.8|) , we find that the transverse components of the vector potential 
read^ 



A^{x) = g / dyii U(?/||,a;) 



— oo 
oo 



d'^d^G{y,~^„x-^)p{^[ 



U-'(2/,„^) 



^^EH^T"' r ^"^U / d'i,d^Giy,^-^,^;x-^,) 

m=l ''-•^ 



^The longitudinal component, A~{x) turns out to be of order e. In any event, it does not 
contribute to the gluon number density, as it does not represent a physical gluonic polarization 
state. 
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\t=2 



n / d%G{y,, - e,,,; x - [[p(fi)p(6) ■ ■ ■ p{i 



(5.13) 



The quantity in the double square brackets appearing in the last hne of Eq. ( |5.13|) is simply 
a multiple nested commutator: 



p(6)p(e2)---p(e« 



p(ei),p(6) ,p(e 



(5.14) 



At this stage there are no more explicit factors of e appearing in the vector potential or 
in the expression for the gluon number density (p.2|) . Thus, the e — limit is trivial to 
perform. 

The expression for the gluon number density involves the partially Fourier-transformed 
quantity A{q^^] x). Since the only dependence on itself appearing in Eq. ( |5.13| ) is as the 
upper limit of the outermost of the ordered integrations, we have 

Ai{q,-x) = ^ E {-igT-' r exp(-zg,||/|, J / d% d^G{y,^ - ^i,,; x - Ci) 



m=l 



' — oo 
m 



II d%G{y,,-^,,;x-^e)] p(6)p(6) ■ ■ ■ p(U • (5.15) 



Finally, we insert the mixed representation of the Greens function presented in Eq. ( ^.6[) . 
The transverse part of the resulting ^ integrations simply Fourier transforms the transverse 
part of the charge densities: 



P(eii;p) = Jd'^e^^-^pii). 



(5.16) 



Hence, the light-cone gauge vector potential becomes 

1 °° roo POO 

A^{q,;x) = ^ E i^9T / d^y,^ / ^-^1, exp(-zg,||/|, J 

IQQw 1 J — oo J—oo 



m=l 



X 



X 



d'^Pl Plj 



^e-*P-^exp 



n 



47r^ 2pi 
d^Pe 1 



e ^^^■"'exp(-p^|y||^-^|| 



2^ 47r2 2pe 
X [[p(^iii; Pi)p(Cii2; P2) ■ ■ ■ P(^ii^; Pn 



(5.17) 



It is useful to have a diagrammatic representation of the contributions to Eq. ( p. 17] ). 
Because it turns out that the longitudinal structure is significantly more complicated than 
the transverse structure, our diagrams are meant as an aid in understanding the longitudinal 
structure. Fig. |] illustrates the mth term of Eq. ( |5.17|) . The vertical line represents the 
range of the y^^ integrals, with each of the vertices (dots) being the value of one of the i/n's. 
Because these integrations are ordered, the dots are not allowed to slide past each other. 
The sources, whose longitudinal coordinates are (^i^^, ^1,2, . . . , ^n^), are denoted by the circled 
crosses. The lines connecting the ?/||'s with the ^n's correspond to the longitudinal factors 
of the Greens functions, exp{—pi\y^^. — ^nj). Finally, the index j labelling the y\\i-^\\i line 
reminds us that it is special: not only is the Greens function associated with this factor 
differentiated (producing the factor of Pij), but there is a factor exp{—iq^^y^^-^) left over from 
the longitudinal Fourier transform which was performed on the vector potential. 
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C. Determining the Gluon Number Density 



We now turn to the computation of the gluon number density. Because of the extended 
longitudinal structure, the calculation is rather lengthy. Here we will outline the path to the 
result with the help of our diagrammatic representation. The mathematical details appear 
in Appendix |^. 

The full calculation of the gluon number density essentially consists of inserting two 
copies of Eq. ( |5.17| ) for the vector potential into the master formula ( p.2| ) for the gluon 
number density, and performing all possible pairwise contractions of the sources in each 
term of the result. We retain only the leading terms in powers of a/R. There are a total 
of four longitudinal integrations per contracted pair of sources: the (unordered) position 
associated with the inversion of the Yang-Mills equation ( [5. 21 ) for each source (represented 
by the circled crosses in Fig. ^ plus the (ordered) integration on associated with the 
transformation to the light-cone gauge (represented by the points on the vertical line). 
A priori, these integrations could produce a factor of (per pair). However, we shall 
now argue that at most they produce a factor of a^R. The key observations to make 
are that both the propagators and the correlation function (pp) (through V) allow for a 
longitudinal separation between the points they connect which is at most of order a. For the 
(pp) correlator, this property follows immediately from the color neutrality condition: the 
confining nature of QCD tells us that T>{x — x') should be negligible when the two points 
being compared are separated by more than the nucleon size a (transversely or longitudinally, 
since we take V to be spherically symmetric). In the case of the propagator, the mixed 
form ( ^.6|) is particularly illuminating: the longitudinal separation of the two points must 
be order 1/p or less (where p is the typical transverse momentum flowing in the propagator) 
to avoid exponential damping of the contribution. However, the (pp) correlation function 
limits the transverse momenta to the region p ^ 1/a, again because of color neutrality. Since 
the four points in question are connected via two propagators and one contraction, three of 
the four integrations are restricted to have range a, while the remaining integration has the 
potential to roam freely over the full range of order R. 

Nevertheless, not all combinations of contractions produces the maximum factor a^R for 
all pairs. Fig. ^ illustrates some of the possible contributions at 8th order in p. From the 
previous paragraph we know that each of these diagrams could, at most, contribute four 
powers of R. However, the contribution in Fig. ^ contains only three powers of R: the 
self-contraction connecting ^||2 with effectively forces y\\^ and j/y^ to be at most a distance 
a apart. But this hems in the point at yyg, preventing it from independently spanning the 
full range R. On the other hand, the self-contraction in Fig. does not restrict the range of 
any additional ?/y 's: it contains the maximum four powers of R. Likewise, the set of mutual 
contractions illustrated in Fig. ^ produce only three powers of R, since the two "crossed" 
contractions cannot slide up and down independently. In contrast. Fig. |^ contributes the 
full four powers of R, since all of the "rungs" may move freely through the full vertical range. 

Thus, the computation in Appendix ^ includes all diagrams which contain only un- 
crossed mutual contractions (like Fig. ^d) or any combination of self-contractions between 
adjacent sources plus uncrossed mutual contractions (like Fig. |^). This produces the leading 
behaviour in the limit a/R -C 1. The final result reads 
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dN 



3ACf'^ — [ (fAe"'-^C(x-,A)S(v^L(A) 

71^ J ^ 



dXp(P'q 



(5.18) 



where 



1 -ip.A p'^V{Xpm,p) 



47r2 



+ {XpTny 



2 ' 



(5.19) 



and 



L(A) 



rd^pV{Q,p) 




J 47r2 p4 





(5.20) 



Despite their superficial appearance, these functions are infrared finite for a spherically 
symmetric correlation function which satisfies the color neutrality condition (|3.7|) . In the 
Abelian limit {a'^A^^^ 0), the gluon number density is simply 



dN 



dXpd^q 



'iACp'^—C^Xp^q). 



Us 1 q'^V{Xpm,q) 



TT" X, 



q^ + {XpTTiy 



2 • 



(5.21) 



The nuclear correction function S encodes the non-Abelian effects and depends on the ge- 
ometry and size of the nucleus: 



S(z) 



In either case, we have 



z 
3 



^ — 1), (cylindrical); 
1 — + e^{z — 1) , (spherical). 



(5.22) 



lim £(z) 



1. 



(5.23) 



Finally, the magnitude of governs the relative importance of the nuclear corrections: 

^ 2^TxalA^!^ k%p^, (cylindrical)- 



(5.24) 



47ri?2 



367rQ;^y4"'^/^AQQj-,, (spherical). 



These corrections are enhanced for very large nuclei. In writing down Eqs. (|5.22|) -( ^.24|) we 
have assumed that the nucleons are uniformly distributed within the volume of the nucleus. 

The result presented in Refs. is recovered in the — limit by using cylindrical 
geometry, since these papers assume that ^ {i.e. the part of the correlation function which 
corresponds to S in the present paper) is a function of the longitudinal coordinate only. This 
is only true for a cylindrical nucleus. Actually, the difference between the two functions in 
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Eq. ( |5.22| ) is very small when the different values of f ^ indicated in Eq. (|5.24|) are taken into 
account. 

Away from Xp = there are two distinct sources of finite Xp corrections: the correlation 
function T>{Xpm,p), and the propagator appearing in the Abelian result. What is perhaps 
surprising is the fact that the function L{A.) turns out to be identical to its 2-D counterpart: 
it depends only on the value of V at Xp = 0. 



D. General Properties of the Gluon Number Density 

Because of the similarity of Eqs. ( |5.18|) - (|5.24|) to the result obtained in Ref. 0, the 
properties of the gluon number density which were described in Sec. IV of that paper continue 
to hold. In particular, the A ^ behaviour of Eq. ( |5.19| ) is unchanged from the behaviour 



of its counterpart in Ref. [y]. Thus, we still have the transverse momentum sum rule 



1^, 



dN 



dXpd^q 



dN 



all 
orders 



dXpd^q 



lowest 
order 



0, (5.25) 



even when 7^ 0. Eq. ( p.25| ) states that the nuclear corrections have no effect on the 
total number of gluons at each value of Xp-. we could have obtained the same number of 
gluons by ignoring the non-linear terms in the light-cone gauge vector potential. What these 
corrections actually do is to move gluons from one value of the transverse momentum to 
another. Thus, the total energy in the gluon field at a given value of Xp is affected by the 
non-Abelian terms. The gluon structure function resolved at the scale is obtained by 
integrating the fully differential number density over transverse momenta satisfying \q\ < Q 
(see Eq. ( |3.4| )). Consequently, the transverse momentum sum rule tells us that for large 
values of Q^, the non-Abelian effects die off, reflecting the expected asymptotic freedom of 
the theory. We should caution, however, that unless Xp is very small, the maximum for 
which our treatment is valid is not very large (see Fig. |l|). Thus, we conclude that at such 
values of Xp the non-Abelian terms are always important at the (smallish) values of for 
which our approximations hold. 

The over-all shape of the fully differential gluon number distribution is insensitive to the 
detailed nucleon structure incorporated in V. Instead, it is fixed only by the confinement 
scale a ~ Aq^q plus the relative importance of the nuclear corrections, governed by f^. 
Recall that according to the discussion of Eqs. ( p.l8D -( ^.24D , the nuclear corrections are 



contained in the function £^(f^L(A)). To understand how £ behaves, we need to know two 
facts about L{A). Firstly, according to Eq. ( |5.2(J| ), L{0) = 0. Therefore, at large q^, the 
all-orders distribution is identical to the lowest-order result. (This is one of the observations 
used in Ref. @ to derive the transverse momentum sum rule.) Secondly, for a wide range of 
physically reasonable choices for T>{x — x'), L{A.) < 0, and decreases as |A| increases [0. 
Thus, we are interested in the behaviour of S for negative values of its argument, which we 
have displayed in Fig. ^. From the figure it is easy to see that the long-distance contributions 
to the integrand of ( ^.181) are damped by the presence of S. This behaviour is consistent 



with the confining nature of QCD, which we incorporated into the form chosen for T) by 
requiring it to obey the color-neutrality condition. Independent of the other details of V, 
we find that 
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lim Xr 



dN 



dXpd^q 



constant. 



all 
orders 



(5.26) 



that is, the distribution saturates as is lowered.^ Furthermore, since increasing the size 
of the nucleus increases f^, which in turn increases the amount of damping provided by £ 
for the same value of A, the gluon density at which saturation occurs decreases as A^f^ 
is increased. Because the large-g part of the distribution does not change, we conclude 
that in order to satisfy the sum rule ( |5.25| ), the number of gluons at intermediate momenta 
must increase. Heuristically, the position space width of the non-Abelian factor £^(f ^L(A)) 
goes like . This width provides a second length scale in addition to the scale Aq^j^ 
characteristic of the lowest-order result. Thus, we might expect that momenta of order 
~ v'^ would play an important role in the resulting all-orders distribution. According to 
Eq. ( ^.24 ), we expect oc A^/^Aq^j-,. This is, in fact, what we observe in our numerical 
calculations: an enhancement in the number of gluons with transverse momenta of order f ^ 
(see Fig. || in Sec. |VT|). Although the idea that a new scale proportional to A^/^h?Q^Qj^ should 
emerge and play an important role for large enough nuclei is not new our 
results lend further support to this concept. 



VI. ILLUSTRATION OF OUR RESULTS 

A. Power-Law Model for V{q) 

In this section we will illustrate the features of the gluon number density described in 
the previous section by choosing a specific form for 'D{q), namely 



1 



[1 



(6.1) 



where uj is an arbitrary positive integer, and = a/y/3uj. This function obviously satisfies 
the color neutrality condition (|3.7|) . In terms of the Kovchegov model |Tl|, the choice made 
in Eq. ( |6.1| ) corresponds to a Yukawa-like distribution of quarks within the nucleon (see 
Sec. [VI Bp . The value of has been chosen so that the root-mean-square radius of the 
nucleon in precisely a. Eq. (|6.1|) is also convenient in that we are able to perform the 
integrals in Eq. ( p.l9| ) and ( |5.2CI| ) analytically (see Appendix 



47r 



[x pVna^jY Kq{x pm/S) 



1 / A y' l-f {uj-j){xpma, 



47r j! V2a, 



1 + {Xpma^y 



[Xpma^ 



(6.2) 



The lowest order result ( 5.21| ) actually vanishes at = when = 0. This may be viewed as 



an accidental cancellation in the integrand of ( 5.18 ) when £ = 1. 
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and 



L(A) = - ^uj 



a,,,/ V2a, 



2 u,-\ 



1 WAV^^AA 



2j j! V2a^/ ■^Vcia;/ 



(6.3) 



In Fig. ^ we have plotted the integrand of the gluon number density ( |5.18| ) (omitting the 
exponential factor) for various values of uj: according to Eq. ( |3.2|) this is just proportional 
to (y4^(g||; a;)74"(— gii; x')). We see from the plots that the non-Abelian corrections become 
more important as lo increases: the range of the (position space) integrand decreases. Thus, 
we would expect to find fewer soft gluons in a model with larger lo. Overall, however, the 
dependence on lo is rather weak. Therefore, we have chosen to present plots only for the 
= 1 case for the rest of this discussion. 
Fig. 1^ contains plots of the fully differential gluon number density as a function of at 
Xp = 0.0 and 0.1.0 Three different values of a'^A^^^ have been used, namely 0.0, 0.5, and 
2.0, corresponding to the Abelian limit, a (roughly) uranium-sized nucleus, and a very large 
(toy) nucleus. In all cases the distributions saturate as — * 0, with the turn-over occurring 
at a few times Aq^^j^. This turn-over is very much like the one which Mueller sees in his 
calculation based upon onium-scattering |]9|,p^. In each plot, the maximum value reached by 



XpdN/dXp(Pq decreases as A^^^ is increased. At large q^, the distributions exhibit the 1/g^ 
fall-off characteristic of individual point charges. Fig. |^ illustrates the same fully differential 
gluon number densities, but multiplied by a factor of q^. These plots are useful because 
the visual area under the curves (using logarithmic horizontal and linear vertical scales) 
faithfully reproduces the result of the integration defining the gluon structure function ( p. 41) : 
what-you-see is what-you-get. From these plots, we see that the very small q^ region makes 
very little contribution to QA^Xp^Q'^). At very large q^ all of the curves converge to the 
same result, as required by the transverse momentum sum rule. We also see a pile-up of 
gluons in the region of a few times Aq^^j^. As the size of the nucleus is increased, this peak 
shifts to larger q^ and increases in size. In Fig. | we track the location of this peak as a 
function of f^. This plot clearly shows that for large-enough nuclei, our heuristic argument 
of Sec. |VD| claiming that Q'^^ak ought to be proportional to f ^ is very nearly correct, with 
a proportionality constant close to unity. When A^/^ is too small, however, this relation 
breaks down as the scale associated with £^(f ^L(A)) no longer dominates the result. 

In Fig. 1^ we plot the gluon structure function per nucleon as a function of Xp for various 
values of Q^. Although we have drawn all of these curves over the entire range from = 
to Xp = 0.25, we remind the reader that the maximum Xp value for which our calculation 
can be trusted decreases as is increased (see Fig. p. Away from the very small Xp 
region, we see that as Xp increases, XpgA{Xp,Q^) decreases, indicating a fall off which is 
more rapid than l/Xp. The degree of dependence on A^^^ goes down as is increased: at 
= 1000Aq(^j3 there is very little nuclear dependence beyond the trivial scaling with the 
number of nucleons. 



^By Xp = 0.0 we really mean some value of Xp in the MV region, Xp A '^/^ /{ma). 
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In Fig. [ly we further explore the nuclear dependence of our result, by plotting the gluon 
structure function per nucleon as a function of a^A^/^ at fixed Xp = 0.1 and several different 
values of Q^. At low Q^, we see a marked departure from the naive expectation: the number 
of gluons is not simply proportional to the number of nucleons. Instead, our gluon structure 
function grows more slowly than A as the number of nucleons is increased. At larger values 
of Q^, however, the nuclear dependence is reduced as the non-Abelian corrections become 
less important. 



B. Connection to the Kovchegov Model 

In this subsection we put our power-law model into the context of the Kovchegov model 



of Ref. |TT[. In this model we imagine a nucleus of radius R which contains A nucleons, 
each of radius a ~ Aq^j-,. Each of the "nucleons" is made up of a quark- antiquark pair. 
Under these assumptions, Ref. [|1T| provides a means for computing the two-point charge 



density correlation function. The correlation function derived in Ref. is applicable 
to the 2-dimensional case. However, since the 2-dimensional form was arrived at from a 3- 
dimensional one by integrating over the longitudinal variables, it is not difficult to modify the 
derivation of Ref. ||TT| to cover the 3-dimensional case we are studying here. Parameterizing 
the correlation function as in Eq. (|3.5|) we find that 



1^/x + a^ Jp(f _ f ) ^ J^ng{x,x') -Ismth{x,x'), (6.4) 



V \ 2 
where 

Isng{x,x') = J d'f\<P{f)\'J d'(m)\'5\x-f-i)5'ix'-f-i) (6.5) 

and 

IsnM^') = I d^rW)? I d'im)\^ I d'i' \i;{^')\^8\x-r~i)8\x'~r-^'). (6.6) 



In Eqs. ( p.5|) and ( |6.6|) the position of the nucleon relative to the center of the nucleus is 
denoted by f while the position of the (anti) quark relative to the center of the nucleon is 
denoted by ^ In addition to making the result fully 3-dimensional, we have allowed for 
the possibility that the quarks and nucleons have some distribution other than uniform: 
the function 10(7^)1^ gives the probability distribution for nucleons within the nucleus while 
|?/'(^)P does the same for the (anti)quarks within the nucleons. Both functions are nor- 
malized to unit total integral. If we compute them exactly from Eqs. ( |6.5|) and ( p.6| ), the 
functions Tsng{x, x') and Xsmth{x, x') do not factorize on the sum and difference variables as 
implied by the left hand side of Eq. ( |6.4| ). However, if we approximate the integrals by as- 
suming that a <^ R (the large nucleus approximation), then they do factorize as advertised, 
and we obtain 

S(^) = Vm)\'; (6.7) 
V{A ) = 5%A ) - I d'i 1^(0 P / d^e \^(i')\^5'ii- f - A ). (6.8) 
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The separation of the right hand side of (|6.4|) into S and V has been uniquely fixed by 
imposing the normahzation for S specified in Eq. ( p.8|) . 

To make the connection with our power-law model, we first Fourier transform Eq. ( |6.8| ) 
to obtain 



1-Viq) 



(fx e''^-- 



[X 



(6.9) 



For a spherically-symmetric nucleus, V^q) is real. Thus, Eq. ( p. 91 ) is telling us that the 
square root of l—T) is given by the Fourier transform of Consequently 



[X 



d?q 



-iq-x 



l-V{q). 



(6.10) 



For the power-law model ( |6.1|) , the integral in Eq. (|6.10|) is easily performed using the 
methods outlined in Appendix M, with the result 



[X 



4(7ra2)3/2 r(../2) V 4a^ 



2 , 2 "-3 



K^i—Jx'^ + x'^ 

2 \a, 



(6.11) 



Eq. ( p.ll| ) describes a Yukawa-like distribution of quarks within the nucleons in the sense 
that the long-distance behaviour of the modified Bessel function is a dying exponential: 



Attx^Wx)? 



exp + x„ 



(6.12) 



At the origin, we have 



constant, cu = 1: 



\x\-.0 



0, 



00 > 2. 



(6.13) 



These distributions are plotted in Fig. 11 for = 1, 2, and 8 



VII. CONCLUSIONS 

We have extended the McLerran-Venugopalan model for the gluon distribution of a very 
large nucleus to larger values of x^. The classical computation contained in this paper, 
like those in Refs. [0-0], is based on the premise that the quantum corrections may be 
ignored since we are working in a regime where <C 1. However, at small x^ there are 
quantum corrections to the distribution functions of the MV model which are proportional 
to as\ia{l/Xp) 0,^. For fixed ~ ^qcd' "^^^ ^e shown that asln(l/x^) is of order 1 
in region C of Fig. and that furthermore, asln(l/a;^) increases as x^ is decreased into 
regions A and B. This suggests that the quantum corrections may be more manageable in 
region C: a leading-log calculation to resum these contributions might suffice. Considerable 
effort P, p4| , p7| -|31[| has already been performed with the goal of dealing with the quantum 



corrections in the very small x^ region where asln(l/a;^) ^ 1 
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The values of we have considered in this paper are in the regime where the gluons 
begin to probe the longitudinal structure of the nucleus. A description of the physics in 
this situation must begin with a fully three-dimensional framework and a source which is 
slightly off of the light cone. We have solved the Yang-Mills equations for such a source in 
the covariant gauge to lowest order in g^, and then transformed that solution for the vector 
potential to the light-cone gauge, where the connection to the gluon number density is to 
be made. The determination of the gluon number density itself relies heavily on the fact 
that we are considering a large nucleus {R ^ a) which consists of color-neutral nucleons. 
We have obtained a relatively compact expression in this limit, which sums the non-Abelian 
effects to all orders in a; 

^1/3 _ Yoi 0, our results match smoothly onto the previous 

treatments Our results for the gluon number density exhibit saturation at small q^: 

instead of diverging as — 0, the distributions approach a finite constant, as illustrated in 
Fig. |. The nuclear corrections induce a pile- up of gluons at ~ w^, where f ^ oc A^/^Aqqj-,. 
In addition, the gluon structure functions which we obtain grow less rapidly than A as the 



number of nucleons is increased, as illustrated in Fig. 10 
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APPENDIX A: EXPLICIT COMPUTATION OF GLUON NUMBER DENSITY 

In this Appendix we will present the details of the derivation leading to the expression 
for the gluon number density presented in Eqs. ( p.l8[ )-( ^.24| ). Our starting point is the 
light-cone gauge expression for the vector potential: 



m=l 



1 POO POO 

~e"""-exp(-p,|!,„-?,j) 



e=2 

X 



[[pten; Pi)ptei2; P2) ■ ■ ■ ptei^; Pm)]] • (Al) 



Simply inserting two copies of Eq. ( |A1|) into the master formula (|3.2| ) leaves us with a 
bewildering array of contractions which must be performed. A more efficient way to proceed 
is to organize the computation in two stages, as was done in Ref. 0. First, we consider all 
of the ways in which pairs of p's within a single (gy ; x) may be contracted. We will see 
that these self-contractions exponentiate provided that a/R <^ 1. Afterwards, we will deal 
with the mutual contractions, where one p comes from each factor of A^{q\'^ \ x). In the end. 
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we retain exactly the same terms as were retained in Ref. |^. However, the justification for 
keeping only these terms is very different from the one in Ref. [^, and relies heavily on the 
large nucleus approximation, a/i? -C 1. 



Since Eq. ([Al|) contains charge densities which have been Fourier-transformed on the 
transverse variables, it is convenient to do the same to the two point correlation function. 
It is straightforward to show that the corresponding transform of Eq. (p.5|) reads 



1. Self Contractions 

In this subsection, we will show that the various non-vanishing self-contractions within a 
single A^{qf x) may be arranged into an exponential factor. We will use the normal-ordered 
notation : ■ ■ ■ : as a bookkeeping device to indicate which factors are not to undergo further 
self-contractions. 

We begin with the observation that the contraction between and p(^|i2;P2) 

vanishes identically: Eq. ( |A2D is symmetric in the color indices, whereas the factors being 
contracted are antisymmetric, appearing in the innermost commutator. 

Next, we show that in the limit a/R <^ 1 {i.e. A^^^ ^ 1), the only non-vanishing 
self-contractions are between adjacent factors of p (such as the self-contraction illustrated in 
Fig. ||b). Suppose we consider a term in which we contract the non-adjacent factors p{C,i\^; Pi) 
with p{i\\pPj), where j > i + 1: an example of such a contraction is illustrated in Fig. |^a. 
The relevant longitudinal factors coming from this type of contribution readQ 

/OO /"OO 

x/^^^iii/^c?^iij exp(-pi|i/||. -ei|.|)exp(-pj|?/ii^. 

where the 0-functions encode the ordering required in the j/n integrations. To see how this 
contribution is subleading when we make the large nucleus approximation, we first introduce 

'^\\i = y\\i-^\\{^ '^\\j = y\\j-^\\j (A4) 

in favor of ^y. and ^y^., followed by 

S|i = -T|u+2/||,-T„^.); All = 2/iu-T|| -2/,i,.+T||^. (A5) 

to replace T/y^ and y^^.. The Jacobians of both transformations are unity. Hence, ( |A3| ) 
becomes 



^The additional factor of exp(— j^) which would also be present \i i = 1 does not affect the 
outcome of our argument. 
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/oo roo roo roo 

rfS,,/ dAj dT,J rfT„^.e(y„„,-S,-iA,-T„je(E„+|A„+T„-y„_^J 
-oo J— oo J ~oo J —oo 

xe(y||^._,-S||+iA,-T||^.)0(S„-iA,|+T,|^-y,i^.^J 
X exp(-p,|T||J-pj|T||^.|) 

X + P,)I^(A,|; - p,)), (A6) 

Recall that the function S{T,) is nonvanishing provided that |ZI | ^ R. Likewise, the function 
P(A) is dominated by the region where |A | < a. Therefore, as far as the transverse 
integrations are concerned, the integrand in ( |A^) is dominated by the region Sy ^ Ay. 
Neglect of Ay relative to Sy in the B-functions will result in errors of order a/R. Furthermore, 
the exponential factor restricts the values of Ty^ and Ty^. for which the integrand is significant 
to Ty^, Ty^. < Pi,Pj- So unless the region where pi or pj is ^ 1 / a is important, we may also 
drop T||^ and Ty^. from the B-function arguments. However, we know that the typical 
momenta associated with 5(Sy;pj + Pj) are \pi + Pj\ ~ whereas the typical momenta 
associated with ^(^Af, — Pj)) are ^\pi — Pj\ ~ 1/a. Together, these constraints imply 
that the main contributions to the integral occur when pi and pj are back-to-back to within 
an amount of order 1/R, and they each posses a magnitude of order 1/a. Thus we conclude 
that Ty^ and Ty^. are indeed of order a, and may be dropped from the B-functions. Making 
these approximations in { \KB^ ) yields 

/oo roo roo roo 

rfSy/ dAj dT,J f/Ty.B(yy,„,-Sy)B(Sy-l/y^^JB(yy^._,-Sy)B(Sy-l/y ) 
-OO J —OO J —oo J —oo 

X exp(-p,|TyJ-p,|Ty^.|)5(Sy;p, + p,)v[A,- \{p,, - p,)). (A7) 

The B-functions in this expression tell us that Sy should lie between y\\^_^ and Vw^j^i on 
one hand, and (simultaneously) lie between y\\j_i and ?/y^^_^ on the other. However, the 
?/y's are ordered, and since we are considering non-adjacent factors, j > i + 1, these to 
ranges do not overlap. Thus, ( |A7| ) vanishes, and we conclude that the contributions from 
contraction between non-adjacent p's are suppressed by one or more powers of a/R relative 
to contributions from contractions between adjacent p's. 

Now that we know which self-contractions may be ignored, let us begin a term-by-term 
examination of the series in Eq. ( [AID • We will denote the mth term in the sum by ^(m) {Qw ', x) . 

The first term, y4|j^^(gy; x), has only a single factor of p. Thus, we trivially obtain 

Al^{qf,x) ^ ■.Al^{qf,x):. (A8) 

Likewise, since the only possible self-contraction which we may consider extracting from 
^^2) {q\i ', x) vanishes, we have simply 

Al^^{q,;x) ^ : Aj2)(gy;a;) : . (A9) 

At third order, in addition to the contribution where we choose to do no contractions, 
we have a term which is generated from contracting p{^\\2-,P2) with p{S,\\^'-,p-i)- The color 
algebra associated with this contraction is straightforward: 
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X( 

X 



The interesting longitudinal factors read 

/oo roo roo roo 

dy\\2 ^^112/ diu^hj\\i-y\\2)^hj\\2-y\\z) 

-00 J ~OQ J —00 J —00 

;exp(-p2|2/|i2 - -Ciial) exp(-p3|yii3 - -^n 
We again make the variable changes indicated in Eqs. ( |A4|) and (|X5|), producing 

/oo /"OO /"OO /"OO 

rfSii / dA, / rfT„, / rfT„3 0(y||^ - S„ - iA„ - T,,,) e(A„ + T,,, - T„3) 
-00 J —CO J —ao J —00 

X exp(-p2|T,|2| -p3|T||3|)5(S||;p2 +P3)I^(A||; ^(^2 -P3' 



(AlO) 



(All) 



(A12) 



We may apply the large nucleus approximation to drop Ay and relative to Sy in the 
first 0-function appearing in ( |A12| ). However, the same arguments which allow us to do so 
also tell us that Ay, T||2, and are all of order a. Hence, the second B-function cannot 
be simplified. Nevertheless, dropping Tyg and Ay from the first 9-function is sufficient to 
allow the integrations on Tyg and Tyg to proceed easily, yielding 

/ dSye(yy,-Sy)5(Sy;p2+P3) / A y P ( A y ; ^ (^2 " P3 ) ) • (A13) 

P2P3 J ~oc Jo ^ ' 

For a spherically symmetric nucleus, D is an even function of Ay. Thus, we conclude that 
the Ay integration simply completes the Fourier transform of 1^/2, with the longitudinal 
momentum evaluated at zero: 



P2P3 



'J —00 



(A14) 



Applying the results in ( |A10[ ) and ( |A14| ), we find that 



^(3) (^11;^) 



^(3)(9ii;*) + 



£pl Plj 



X 



47r2 2pi 



exp(-igy?/y^ 



(A15) 



In (|A15|) we have introduced the function 



r(.,;.,.')^£dE,/^/ 



(Pk r (Pk' e-^'^-^ e"^*''-^' 



47r2 J 477^ fc2 



fc'2 



which will prove to be useful as we proceed with the calculation. 

At fourth order, there are two different non-vanishing contractions. Their contributions 
differ only in the range of the Sy integration, and combine neatly to produce 
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^(4)(?lh ^) • 
I ^ ^2 



2\2 



d pi f d^p2 pijC e 



47r2 y Att^ 2pi 2p2 
X / ci^^ii exp(-pi|t/||^-^||^| -P2|yii2-Cii2l) 



X : 



Finally, at fifth order, in addition to the three different ways to perform a single contraction, 
we encounter a contribution containing two contractions. The result of a straightforward 
computation is 



^(5) ; ^) ^ ■ ^(5) iiw ; ^) ■ 



(5)^ 



2\3 



d^pi f d^p2 f d'^ps pije-'P^ '' e-'P2 =" e-'P^ '^ 
4^ 



f g P2 f 

J 47r2 J 



47r3 2pi 2p2 2p3 

3 



/oo roc / ^ 

Alu exp(-ig|||/|| J / rf^^ii exp - ^p^li/ii^ - 
-oo J —oo \ 

X : [[p(^iii;Pi)p(^ii2;P2)p(^ii3;P3)]] : {-|A^c^^X'r(y||,;a;,a;)} 



d'^Pi pij 



/oo 
dy \ii exp(-ig||y|| J 
-oo 



X 



X 



47r2 2pi 

/ c^Cii 1 exp (-pi 1 2/11 1 - Cii 1 1 ) : [UCii i ; Pi ) 

J — oo ^ / LL 



2! 



(A18) 



At this stage we can see the pattern which is emerging. When we choose to do no 
contractions, we get back the series for A^{qf,x), but normal-ordered. Starting at third 
order, we have the option of doing at least one contraction. Choosing to do exactly one 
contraction at each order produces a series which is nearly the one for A^{q^^; x), but with 
an extra factor 



lN,g'K'r{y,^;x,x) 



(A19) 



inserted into the integrand of each term. At fifth order, we may elect to do at least two 
contractions. Doing exactly two contractions in each term again nearly reproduces the series 
for A^{q^^; x), but this time with an extra factor 



^{-^N,g'K'T{y,,;x,x)y 



(A20) 



in the integrand. In like manner, all of the terms in which we do exactly j contractions sum 
up to (nearly) produce the series for A^{qf, x), but with the extra factor 
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j[-\N,g'K'T{y,^-x,x)Y. 



(A21) 



Thus, we conclude that systematically accounting for all possible self-contractions results in 



2\m 



d^Pl Plj 

47r2 2pi 



n 



47r2 2pe 



X 



X 
X 



/oo roo / 

-oo J —etc \ 

p(^iii;Pi)p(^ii2;P2)---p(^ii^;p„ 

exp{-lN,g^K^T{y,,;x,x)}. 



(A22) 



2. Mutual Contractions 

We now insert the required two copies of Eq. ( |A22| ) into Eq. ( |3.2| ), the formula for the 
gluon number density. Because all of the self contractions have already been accounted for, 
we may only multiply terms which contain the same number of p's, leading to a single sum 
(rather than a double sum). Thus, we obtain 



dN 



1 1 



n 



-i//./.v.-<-.i:(-/)™/^/ 



rfVl [ d^P'l Pi ■ P'l ^-ip^.a^-ip'^.a^' 



m=l 



X 



2^ 47r2 y 47r2 Ap^p'^ 

/OO poo \ 1 



1 AT ^4 r^3 



xexp -iiVe/i^^ r(?/||,; 03, a;) + T{y;^; x', x'] 



(A23) 



Using an argument exactly analogous to the one in Eqs. (|A3|)-( [A7D it can be shown that 
the only non- vanishing contribution to leading order in powers of a/R is obtained by per- 
forming "corresponding" contractions (like the contractions in Fig. ^), that is, p{^\\pPj) 
with p(^||'.;p' ) for all j. "Crossed" contractions (like Fig. ||c) are suppressed by one or more 
factors of a/R. 

The color algebra associated with the corresponding contractions involves the expression 



T = Tr(^r" 



(A24) 



The required sums are easily evaluated by beginning with the innermost commutator: 



a/3 
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(A25) 



The result of inserting ( |A25|) into ( |A24|) has the same structure as we started with, but with 
one less commutator. Repeating this process until just two color matrices remain in each 
trace and doing those traces yields 



m— 1 



T=\{Nl-l){-N^) 
Now let us consider the longitudinal integrations. The relevant factors are 



(A26) 



X 



dry, / dry; / d^i, / d^^; exp -^^,,(^1,^ -^„;) 

oo J—oo J—oo J ~oo 

m 

n Qiy^.^-y,,) 0(yi,;--i-yii;) s^IH.^HQ-.p^+p'^ ^^(eii,-^,,;.; \{p,-p-)). (A27) 



where we have defined y^^ = y^^ = oo for convenience. We perform variable changes which 
are completely analogous to those in Eqs. ( |A4D and ([ASD and once again apply the large 
nucleus {R ^ a) approximation: 



(i"E|| / rf'"A|| / c/™T|| / rf'^r; exp -2g,|(A,|,+T||,-T||;) 

) J—oo J—oo J—oo 

X expj - f;(p,|T,|,| + p^|T,|;|) I exp\-lN,g'K' [r(S„,; x, x) + r(S|,,; x' , x'] 



X n e(s,|^_^-S|,^.) 5(S|,^.; D(A,|^- 

The necessary Ty and T'| integrations are all easily performed using 

2p 



(A28) 



rfT|| exp(-ig||T|| -p|T|| 
The 6 integrals simply finish Fourier-transforming P on the longitudinal coordinate: 
|_"rfA|, exp(-2g||Ai|)P(A|,;i(p-p')) = P(g||, " P')) • 

Applying these considerations to ( |A28| ) produces 
2pi 2p[ 



(A29) 



(A30) 



exp|-iAr,(7^i^3 [r(S,|,; x, x) + r(S„,; a;', a^') 
4 ~ 



(A31) 



Notice that the S|| integrations have inherited the ordering associated with the gauge trans- 
formation. When we insert ( |A31| ) back into Eq. ( |A23| ) , the portion of the resulting ordered 
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integrand involving S||2, . . . , Sy^ is symmetric, allowing us to do the sum on m to obtain an 
exponential. Including the color factor contained in Eq. (A2C) we arrive at 

dN 



dq^icpq 



7r2 J J 



dEn 



d^p 



47r2 J 47r2 (pj + qfMp? + q''] 



S{J:,^;pi+p[)v[q„l{p,-p[] 



X exp N.g'^K' [r(S„,; x, x') - ir(S„,; x, x) - |r(S„,; x\ x') 



(A32) 



To proceed further requires us to apply the consequences of the large nucleus approxi- 
mation to the transverse coordinates. To see how this works, let us examine the function F 
a bit more closely. From Eq. ( [Al(j| ) we may write 



d^p e-*P-(^-^')r'(0,p) 

{p + W{p-\q) 



(A33) 



where we have changed variables to q = k+k' and p = \{k~k'). Recall from the discussion 
in the paragraph following Eq. (A6) that the values of the momenta associated with S 
are g ~ l/i? whereas those associated with V are p ~ 1/a. This suggests that we may 
neglect q in the two denominators of Eq. (|A33|) , the error being suppressed by a factor of 
a/R. However, we must be careful. The combination appearing in the square brackets of 
Eq. (|A32|) can be shown to be infrared finite provided that T) is rotationally invariant and 
satisfies the color neutrality condition. This is true to all orders in a/R. When dropping 
terms which are higher order in a/i?, we should avoid introducing an infrared divergence, 
since none was present in the original expression. Therefore we write 



r iy II ^ X ^ X 



d'q 



47r2 



rd'p V{0,p) 






J 47r2 p4 







that is, when we drop q from the denominators we should simultaneously adjust the expo- 
nential multiplying S to be identical in all three terms. The advantage of the form contained 
in ( [A34D is the decoupling of the two momentum integrations. The integral on q just converts 
q) back to a purely position-space quantity. The p integral defines the function 



L{x) 



fd^pV{0,p) 


^-ip-x _ ^ 


J 47r2 p4 





(A35) 



Thus, 



F(S||-j^; X, X 



J — CO ^ 



X 



+ 0[-]. (A36) 



Treating the pi and p'^ integrals of Eq. ( [A32| ) in a similar manner and applying Eq. ( |A36|) 
yields the expression 
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dN . 2a, 1 



X exp|ArcC/^ir3L(a;-a;') ^ + a;')) (A37) 



where we have introduced the quantity 



1 fd'p^^,^.^ p'V{q,, p) 

4:71^ (p2 + g2 



/:(g„;x) = ;/^e---^:^^. (A38) 



Finally, we apply the chain rule to do the integral over Sy^^, and switch to sum and difference 
variables for the x and x' integrations: 

dN _o,„o 2a, 1 



K\N^ - 1) ^ - / rf2Ae'9-^£(g„; A) 

Vr"' On J 



, ^ explg^N.K'^LiA) JZo d^\ S)| - 1 



3. Geometric Dependence 

In order to perform the S integration appearing in Eq. (|A39|) , it is necessary to specify 
the geometry of the nucleus. We will consider two cases, cylindrical and spherical. 
A cylindrical nucleus is described by the function 

= e(i?2_s2)e((M)2-s2), (A40) 



where R is the radius of the cylinder and h is its height. Actually, the height will drop out 
of the final result, since ( |A39|) depends on 



/oo 
-oo 

The 5] integral which results from inserting (|A41|) into ( |A39|) is trivial, producing 
dN 2a 1 /• o ■ . exp\\3Agy27iR^]L(A)] -1 



dq^^d^ " 7r2 ' [3Agy27rR^]L{A) 

which is equivalent to the portions of Eqs. ( p.l8| )-( [5.24| ) pertaining to cylindrical geometry. 
Turning to the more realistic case of a spherical nucleus, we have 



/oo /"OO 
-oo J —oo 



2Vi?2-S2 0(i?2_5]2), (A43) 
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so that Eq. ( |A39| ) becomes 
'^^ K\Nl - 1) ^ 1 / rf^Ae^^-A^^ ^) 



g^NcK^L{A) 



Thus, the S integral hinges upon the form 

rR 



This integral is easily performed by the change of variables 

s = nVR"^ - S2; sds = -n^j: ds. 

Then 



fi2 



Applying (|A47|) to (|A44|) leads to the result 

dN _ „ 2a, 1 



dq^^d^q 



3ACf'^- [ d^Ae^'^-'^Ciq.iA) 
qj 



(A44) 



(A45) 



(A46) 



(A47) 



l-l[v'L{A)r + 



v^L(A) - 1 



exp 



v'L{A) , (A48) 



where f ^ = 9y45('^/47ri?^ 



APPENDIX B: CALCULATIONAL DETAILS FOR THE POWER-LAW MODEL 

1. Useful Integrals 

All of the integrals required to compute the functions C{qf,x) and L{x) which appear 
in the integrand for the gluon number density for the power-law (Yukawa-like) model intro- 



duced in Sec. VI may be derived from the forms 



-iqx 



d^q 

4^ (q2 + g2) [1 + ^2(^2 + ^2)].^ 



(Bl) 



and 



I, 



2oj 



-iqx 



<Pq 

4^ (g2 + g2)2[l + 52(^2 + ^2)].- 



(B2) 
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Since the uj ^ limits of ( |B1|) and ( P2|) are smooth, we may simply set = to obtain the 
necessary single-denominator integrals. Because the procedure for performing both integrals 
is essentially the same, we will describe the computation for and simply quote the result 
for Ji^. 

The computation of l2uj begins by introducing a single Feynman parameter to combine 
the two denominators: 



Z 



2w 



Uj{ijJ + l) 



dz z^-\l-z) 



-iqx 



(B3) 



In order to deal with the q integration, we introduce a Schwinger parameter to promote the 
denominator into the exponential: 



1 



2ui 



(cu - 1)! Jo 



dX X 



OJ+l 



dzz'^-'il-z) 



47r2 



g-iq-cCg-A(q^+g|+2/fe^) 



(B4) 



The q integration is now Gaussian, and may be performed in the usual fashion, with the 
result 



1 



1 



2w 



47r62^ (uj 



—7 / dX A'^ exp 
I)' Jo 



2u 



1 r°° 
— / dX 



4tt Jo 
1 

47r 



1 X2 




A T " 




g 






exp 



f^dzz^-^l-z) e-'^/'' 
Jo 
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1 dX 



e=o 



■aT-^^" 









)1 


exp 


-aT-^^". 



(B6) 



It is convenient at this stage to introduce the dimensionless integration variable ^ = X/b"^. 
Doing so and performing a bit of algebra we arrive at 



2lu 



b^ 






47r2 


/ exp 

Jo 






dC 




47r2 


UJ / — exp 
Jo 4 


62 

^ 47r2 


^ d J 

1=0 ^- 


roo 
D 



1 / X 



^ V26 



rfer-'exp^--(^-) -e i + 



(B7) 



The ^ integrals may be performed to produce modified Bessel functions, as seen from 
Eq. (3.471.9) of Ref. 



(B8) 



which is valid for all values of ;U, provided that (3 and 7 are positive. Thus, we arrive at 



2'. 



2w 



62 r X 1 

2^ 1 26 H 



i^'i(xgii) - ujKo{xq\\ 
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The analogous treatment of (|BT| ) yields 



27T 



A"o(xg|| 



E 



(BIO) 



2. Computation of £(g| ;a;) and L{x) 



Eq. ( |5.19D defines the function which governs the Abelian portion of the integrand for 
the gluon number density. For the power law model, it is helpful to rewrite the numerator 



using = + - gif: 



£(gii;a:;) 



4vr2 



-tpx 



2\2 



1 - 



(Bll) 



In terms of the integrals ( |B1| ) and ( |B2| ) introduced in the first part of this appendix, we have 
simply 



C{qf, x) = iJio - ig^Jao - + ig'Xa^- (B12) 
A straightforward substitution of the results contained in Eqs. ( p9|) and ( piO|) leads to 



UJ 



+ 



1 1 



E 



47r j\ \2a, 



X y 1 + {uJ-J){a^q^^) 



1 + {acuqii 



X 



1 + {au^quy 



(B13) 



The determination of L{x) from Eq. ( |5.2(]| ) is a bit more involved. The difficulty lies in 
the fact that it is not possible to integrate ( 5.20 ) term-by-term, as the individual bits are 
infrared divergent. To work around this difficulty, let us define the auxiliary function 



Lsiqii^x) = 



-ipx 



47r2 (p2 + g2-)2 

Then, the integral we seek may be determined from the relation 



Lix) 



lim 



L^{qi,x) - LsiqfO) 



(B14) 



(B15) 



The computation of -^^3(511; x) is straightforward: inserting the power-law form of V given in 
Eq. (|6.1|) we see that 



LaiQw^x) = J20 -Xa^. 



(B16) 
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Thus, the apphcation of Eq. ( pQ]) gives 

^2 



■ UJ 



27r 

27r ^ 



X 



1 + {acoQwy - KQ{xq^^ 



-112 ( X 



1 + {a^q\\ ) KA — Jl + (a^^gii ) 



(B17) 



To determine L3(g||;0) from (pi7| ) we require the following forms of the modified Bessel 
functions for small values of the argument: 



-In 



^ ^-0 r(/i) /2N^ 



- 7i?, = 0; 



where 7^; is Euler's constant. Consequently, 



^3(gii;0) = ^a;ln 



1 + (a^gii' 



1 



47r,tt ^ [l + (a..g||)2]^ 



Subtracting ( |B19| ) from ([B17|) and taking the gy — >^ limit gives the final result 

^2 



iTn — +ln 



X 



,2 <^-l 



27r 



2a^ 
1 / X 



2£ £! \2a. 



+ 'lE 



(B18) 



(B19) 



(B20) 
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FIGURES 





FIG. 1. Approximate region where the density of color charge is large (and hence weak), 
and sufficient color charge is being probed to justify a classical approximation to the quantum 
average in Eq. ( p. 21 ). (A) Small-x^ region of the original MV model (B) Additional allowed 

region at small Xp when the effects of confinement are included Q. (C) Extension to larger Xp 
discussed in this paper. The upper plot is for a toy nucleus with A^/^ = 250, while the lower plot 
is for uranium-238. 
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FIG. 2. Diagrammatic representation of longitudinal structure of the mth-order term in the 
expression for the light-cone gauge vector potential, Eq. ( [5.17| ). The circled crosses denote the 
positions at which the sources are being evaluated. The dots represent the ordered integrations 
coning from the gauge transformation into the light-cone gauge. 
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FIG. 3. Some of the contributions to the gluon number density at 8th order in p. (a) Diagram 
containing a non-adjacent self-contraction, (b) Diagram containing an adjacent self-contraction, 
(c) Diagram containing a pair of crossed mutual contractions, (d) Diagram containing only corre- 
sponding mutual contractions. Diagrams (a) and (c) are suppressed by a power oi a/R relative to 
diagrams (b) and (d). 
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FIG. 4. Plot of the nuclear correction function £{z) for negative values of z and a spherical 
nucleus. 
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FIG. 5. Position space correlation functions used to determine the fully differential gluon 
number density from Eq. (5.18) evaluated using the power-law model for 'D[q) given in Eq. ( |6.1D . 
The three curves compare the results using w = 1, 2, and 8 at fixed Xp =0.1 in the Abelian limit 
(a2^V3 = 0) and for uranium {a^A^^^ = 0.5). 
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FIG. 6. The fully differential gluon number density Eq. ( |5.18| ) evaluated in the power-law model 
with uj = 1. The three curves on each plot represent the Abelian limit {a^A^^^ = 0), uranium 
{a^A^/^ = 0.5), and a large toy nucleus with A ~ 15000 (q^^^/^ = 2.0). The upper plot illustrates 
the results for Xp — > 0, whereas the lower plot has been drawn for Xp = 0.1. 
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FIG. 7. The fully differential gluon number density Eq. (|5.1 multiplied by and evaluated 
in the a; = 1 power-law model. These plots accurately reflect the relative contributions to the 
gluon structure function 17^1 (x^, Q^) coming from each value of . The three curves on each plot 
represent the Abelian limit {p?^A^I'^ = 0), uranium (a^A^/^ = 0.5), and a large toy nucleus with 
A ~ 15000 (a2^^/3 = 2.0). The upper plot illustrates the results for Xp 0, whereas the lower 
plot has been drawn for =0.1. 
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FIG. 8. Value of the momentum-squared at the peak of the (^dN jdXj^Sq distribution as a 
function of the scale set by the non-Abelian corrections, oc J^l^b^Q^i- These curves have been 
generated within the u = 1 power-law model. 
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FIG. 9. Gluon distribution function gAix^^Q"^) in the uj = 1 power-law model plotted versus 
Xp for Q"^ = Aqqj), IOAqqjj, IOOAqqj-,, and IOOOAqqj-,. The three curves at each value are for 
the Abelian hmit {a^A^/^ = 0), uranium {a^A^/^ = 0.5), and a large toy nucleus with A ~ 15000 
(a2^V3 = 2.0). 



43 



2.000 1 I I I I I I I I 
= lOOOA^ 

1.000 — ' 




0.010— "^-^ 

: xp = 0.1 : 

Q QQ5 I I I I I I I I I I I I I I I I I I I 

0.0 0.5 1.0 1.5 2.0 

FIG. 10. Nuclear dependence of the gluon distribution function gAiXpjQ"^) in the oj = 1 
power-law model for fixed = 0.1 and = ^qcd^ ^.SAqq^, IOAqqj-,, 25Aqqj^, IOOAqqj-,, 
and 1000Aq(-.j). These functions grow more slowly than A as the number of nucleons is increased. 



44 




FIG. 11. Quark probability distributions as determined in the context of the Kovchegov model 
|11| for the power-law correlation function (|6.1| ). The three curves compare the results for to = 1, 
2, and 8. 
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